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Abstract. We present a simple introduction to the techniques of effective field theory (EFT) and 
their application to QCD. For problems with more than one energy scale, the EFT approach is a 
useful alternative to more traditional model-building strategies. The most relevant such problem 
for this discussion is that of making contact between QCD and the hadronic phase of matter. As 
a simple example, an EFT calculation of the bound states of hydrogen within QED is sketched. A 
more significant demonstration of the power of EFTs, the construction of the chiral Lagrangian and 
chiral perturbation theory, is also included. The results provide us with the road map to a complete 
QCD-based theory of nuclear matter at nonzero temperatures and densities, a vital component to a 
quantitative understanding of the phase transition from hadron gas to quark-gluon plasma. 



INTRODUCTION 

Effective field theory (EFT) is a theoretical prescription for constructing theories span- 
ning multiple energy scales. The physics of a system may appear radically different at 
different energy scales, due to low-energy restrictions on available degrees of freedom 
and symmetries. When trying to construct a theory which spans energy scales, tradi- 
tional methods of physics can therefore be difficult to apply. Rather than stumbling on 
this obstacle, however, EFT provides a method to use the physical difference between 
energy regimes to advantage. We present an introduction to the basic ideas of this useful 
technique, from the perspective of its application to nuclear physics. The emphasis is 
pedagogical. Comprehensive and up-to-date reviews can be found in Ref. [1]. For more 
extensive lectures with a similar perspective, the reader should refer to Ref. [2]. 

In a typical experiment designed to produce a quark-gluon plasma (QGP), two heavy 
nuclei collide at relativistic speeds. The experimental signal that results is necessarily 
dependent not just on the nature of the plasma, but on the nature of the hadronic phase as 
well as the physics of the phase transformation. Trying to understand such an experiment 
with only an understanding of the plasma phase, therefore, is as doomed to failure as 
attempting to understand the latent heat of a water-to- steam transformation using only 
the kinetic theory of gases. 

To remedy this problem, we must seek a true quantitative understanding of the physics 
of the hadronic phase. The theory we seek should furthermore be based on the micro- 
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scopic theory of QCD. The development of theoretically-sound hadron interactions has 
been the major problem of nuclear physics since its inception. This theory must be valid 
over the significant temperatures and densities of the phase transformations typical in 
experiment. However, even a more limited approach that is restricted to the hadronic 
phase can be useful. One would like, for example, to estimate the position of the QGP 
transition line. Moreover, it is believed that the phase diagram of QCD presents other 
interesting features such a liquid-gas transition at lower temperatures and at densities 
close to that of equilibrium for cold nuclear matter. 

At zero baryon density, an expansion in temperature T of the free energy density of 
a pion gas can be derived in EFT and reveals the presence of a transition temperature 
Tc at the point where the expansion diverges [3]. A more precise estimate can then be 
obtained by considering the effects of higher-mass mesons [4], in good agreement with 
lattice results [5]. 

At non-zero baryon density p, on the other hand, the situation is much more com- 
plicated. Lattice QCD simulations are hampered by the infamous sign problem. EFT 
might therefore be the only way to study the problem, but it requires a conceptual leap 
from the work with mesons alone, because nucleons bind. Perturbative expansions for 
the energy density of dilute boson and fermion gases are well known: see, e.g., Ref. [6] 
for r = 0. These are expansions in pa^, where a is the scattering length, which is essen- 
tially the scattering amplitude at zero energy. The problem is that the nucleon-nucleon 
(A^AO scattering length is large, so that the corresponding expansion fails way before the 
QGP phase transition. Indeed, the interesting physics of bound states such as nuclei is 
associated with p^/^ ^ I /a. The failure of a perturbative expansion in density implies 
that we must effectively resum such an expansion. 

The reader may be somewhat surprised by the apparent failure in this regime of the 
perturbative approach, upon which so much of particle physics is based. This circum- 
stance may be attributed to the complexities of the nonzero temperature and density 
requirements, as well as the striking disparities of the two energy scales we would like 
to connect. On the one hand, the consensus of the majority of the nuclear-physics com- 
munity holds that in nuclei 

• nucleons are non-relativistic; 

• they interact via essentially two-body forces, with smaller contributions from 
many-body forces; 

• the two-nucleon interaction generally possesses a high degree of isospin symmetry; 

• external probes usually interact with mainly one nucleon at a time. 

By contrast, in QCD 

• the u and d quarks are relativistic; 

• the interaction is manifestly multi-body, involving exchange of multiple gluons; 

• there is no obvious isospin symmetry; 

• external probes can, and often do, interact with many quarks at once. 

It should not be surprising, then, that some new ideas are required to merge these two 
extraordinarily different bodies of theory. Of course, we expect that QCD encompasses 



the physics of hadronic interactions. The root of the problem must therefore he in the 
difference of energy scales. 

In fact, constructing a QCD-based theory of the hadronic phase is a problem which 
involves three separate energy scales spanning three orders of magnitude. The first and 
most obvious to a reader well versed in high-energy physics is the typical energy scale 
ofQCD, 

Mqc^ ~ 1 GeV. (1) 

The masses of all hadrons except the pion fall within this scale ^, and the scale of chiral 
symmetry breaking is thought to be = Anf^^, where f„ ~ 93 MeV is the pion decay 
constant. The second scale, 

M„„c~100MeV, (2) 

represents the typical momentum of nucleons in a nucleus: the inverse root-mean- square 
charge radius of light nuclei, or the Fermi momentum of equilibrium nuclear matter. It 
contains also the pion decay constant itself, the mass difference between the delta isobar 
and the nucleon, and the mass of the pion. The final energy scale is the typical energy 
scale of a nucleus, 

"""^ ~ 10 MeV. (3) 



Mqcd 



The binding energy per nucleon of a nucleus is typically a few MeV. For example, the 
binding energy of ^He is 28.296 MeV, and the binding energy per nucleon of infinite 
nuclear matter is 16 MeV. 

Since the goal is to construct a theory valid over an energy range of three orders of 
magnitude, and spanning regimes in which the physics is quite disparate, the problem 
is somewhat daunting. We might hope to seek inspiration from an analogous problem 
in atomic physics, the simple case of a hydrogen atom in the ground state. A qualitative 
argument suffices for this discussion. We know from experiment that the mass of an 
electron is 

me~0.511MeV. (4) 

The Hamiltonian is, to a good approximation, H = p^/lnig — a/r, where a — e^/4K ~ 
1/137 is the fine-structure constant. From the uncertainty principle we know that the 
electron momentum is p ~ l/R, if the size of the atom is R. Minimizing E{R) = 
1 /ImgR^ -a/R yields R=l /amg, or 

/7 ~ amg = 3.6 keV. (5) 

Finally, the binding energy B — —E{1/ amg) is 

5~ ^aV = 13.6eV. (6) 

This problem, like ours, possesses three distinct energy scales, on the order of m^, p, and 
B. It is easy to see from the above, however, that in the atomic problem all three scales 
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are coupled by the small fine-structure constant a. It is not so clear what the analogous 
coupling constant should be in the nuclear problem. In fact, the lack of a clear small 
coupling constant is one of the main difficulties of nuclear physics. Even though our 
problem involves QCD, the strong coupling strength is not useful to us, for at the 
highest energy scale in the problem ai(lGeV) w 1 and perturbative QCD, valid at larger 
energies, breaks down. 

If all obvious coupling parameters are of order one in the regime of our problem, we 
have no choice but to seek an alternative formulation. What is required is a theoretical 
framework that has the ability to sensibly and efficiently deal with multiple energy 
scales, as the perturbative approach does, but one that also does not rely ab initio on the 
existence of a small coupling constant. The framework which possesses both of these 
characteristics is EFT, as we shall see. 



EFT is a technique for developing theories of problems with multiple energy scales. It 
is applicable in situations where we wish to understand the physics at some low-energy 
scale as the limiting case of a more general problem whose full features are apparent 
only at some higher energy. 

For simplicity, consider a generic problem with two energy scales. The complete 
physics of the problem throughout the full spectrum of energies can be said to be 
described by some Lagrangian density £ (^) in terms of some degrees of freedom (p. At 
a characteristic energy scale E^^^^^, the full properties of L are vital to an understanding 
of the physics. In general, we may or may not know L. Even if we know L, we may or 
may not be able to solve it for the dynamics of the underlying theory. In our problem 
of the hadronic phase, for example, L is the QCD Lagrangian, and the underlying, full 
physics would be QCD. Processes explicitly involving quarks and gluons are seen to be 
important at energies of order E^^^^ ^ ^qcd- significantly lower energies, however, 
these processes "freeze out" and only hadronic degrees of freedom are available to the 
system. In the EFT, we therefore set some energy cutoff A which divides energies of 
order E^^^^^ from energies at which some of the freedoms of the full Lagrangian can be 
neglected. 

The 5-matrix contains all of the information relating initial states to final states in a 
many -body problem. The elements of the matrix can be calculated from the path integral 



where D is the number of spacetime dimensions in the problem. This expression is based 
in the full physics of the elementary, microscopic theory contained in L. 

We can begin to exploit the differences in physics between the two energy scales by 
using A to relabel the fields according to momentum. 
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and integrating over the degrees of freedom in 9^. This gives us 
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where L^^f is defined by 




(10) 



Si^j-j- is obviously a function only of 9^, since we have integrated over <Pjj. We can create 
a series representation for Eq. (10), 



In general, Eq. (11) is a pure mathematical statement without any real physical content, 
but it becomes more meaningful when we make some statements about the properties of 
the operators O^- and coefficients g^-. 

The operators 0^ can in general be quite complicated. We can see, however, that 
they must possess two important physical properties. The first is that, although they 
contain an arbitrary number of derivatives, the O-'s are in fact local operators in the 
sense that they involve fields at the same spacetime point. Nevertheless, their necessary 
dependence on derivatives follows from the fact that they depend only on fields with 
momentum components below A. As such, the uncertainty principle tells us that the 
operators can probe length scales only down to a minimum distance ~ 1/A, and so high- 
order derivatives will be necessary to perform the averaging over smaller length scales. 
The idea is the same as that of a multipole expansion in elecfrodynamics. The second 
observation we can make regarding the O^ 's is that, for an appropriate decomposition (8), 
they must possess all of the symmetries and transformation properties of the underlying 
high-energy theory. Even if a particular symmetry is broken, it will manifest itself in the 
same way in the effective Lagrangian. 

The coefficients can likewise be put in perspective with some physical arguments. 
It is clear that they reflect the "freeze out" of the high-energy degrees of freedom below 
the cutoff A. Therefore, it is true that their particular form is dependent on the nature of 
the underlying theory and the structure of L {q>j^, q>fj) . Furthermore, they are obviously 
a function of the cutoff parameter A. Note, however, that the specific dependence of 
the g- on A is constrained by the familiar principle of renormalization-group invariance: 
physical observables cannot depend on the cutoff because the choice of the latter is 
arbitrary. 

The physical motivation for EFT should now be apparent. It is an extremely common 

circumstance in physics that we wish to treat some low-energy limit of a fundamentally 
higher-energy problem. In such a case, we expect that the low-energy physics will prove 
to be some reflection of the full problem, but with the phase space restricted such 
that high-energy excitations or degrees of freedom are not accessible. Indeed, this is 
arguably the case in any physical problem short of a theory of everything: there exists 
some smaller-length and higher-energy scale at which new physics occurs, but whose 
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details are unimportant to a description of a particular system. Were this not the case, 
physics as we know it would be impossible. The power of EFT, then, is that not only 
does it provide a framework to unify a problem with the energy scale above it, but it 
turns the energy difference from an obstacle into part of the solution: it recognizes this 
principle of limiting the system's freedom at low energy and incorporates it directly into 
the theoretical strategy. 

The EFT is useful even when we know and can solve the underlying theory at higher 
energies. In this case, it is possible to derive the coefficients from our knowledge 
of the underlying Lagrangian ^{(p). The advantage lies in the reorganization of the 
dynamics in the EFT, as it is often the case that the effective degrees of freedom are 
"collective" excitations of the underlying degrees of freedom. This is the case whenever 
there is spontaneous breaking of a continuous symmetry and Goldstone bosons appear. 
An example is magnons, which are a better way to describe low-energy excitations in 
a system with spontaneous magnetization than individual magnetic moments. Another 
example are pions in QCD. 

But EFT is a necessity when either the underlying theory is not known or, as in the 
case of QCD, it is not currently solvable. In this case, we can rely on Weinberg's "the- 
orem": "if one writes down the most general possible Lagrangian, including all terms 
consistent with assumed symmetry principles, and then calculates ^-matrix elements 
with this Lagrangian to any order in perturbation theory, the result will simply be the 
most general possible S'-matrix consistent with analyticity, perturbative unitarity, cluster 
decomposition and the assumed symmetry principles" [7]. There is no known general 
proof of this "theorem", although it has been proven for a scalar field with Z2 symmetry 
in Euclidean space [8]. Nevertheless, it is certainly plausible, as it is really stating that 
the most general quantum field theory is a direct consequence only of analyticity, unitar- 
ity, cluster decomposition, and symmetry. This, combined with the fact that there are no 
known counterexamples, suggests that we would do well to consider the ramifications 
of the "theorem" regardless of its current lack of rigor. 

Weinberg's "theorem" allows us to formulate a plan to solve problems with an un- 
known or insoluble Lagrangian by using EFT. 

1 . Identify the relevant degrees of freedom and symmetries of the problem. 

2. Construct the most general Lagrangian consistent with these limitations. 

3. Do standard quantum field theory with this Lagrangian. 

"Standard quantum field theory" consists of computing the contributions from all dia- 
grams with momenta Q < A, and then renormalizing the result to relate the coefficients 

to the physical observables of the problem. After renormalization, observables should 
be independent of A. According to Weinberg's "theorem", these steps give the dynam- 
ics of the general system, and the process of building an EFT will have introduced no 
spurious information. 

An important issue is that of ordering the infinite number of contributions to any 
observable, known as "power counting". The minimal assumption is that, barring a sup- 
pression by some symmetry, the coefficients in the expansion will be roughly of 0(1) 
once expressed in units of E^^^^^ according to dimensional analysis. This is called the 
assumption of "naturalness", and its validity is a consequence of our having chosen an 



appropriate cutoff A. The perturbative result of this process is a controlled expansion in 
energy, E/E^^^^. However, since observables are cutoff independent, the same expan- 
sion is valid for any cutoff. As we vary the cutoff, strength moves from one contribution 
to another that appears at the same order in E/E^^^^^. So, regardless of the cutoff, to 
any given order in E/E^^^^^ only a finite number of ^-'s need to be considered. One can 
use a finite number of experimental data as input to determine these gfs, and then use 
the known ^-'s to predict everything else, with an error given by the estimated size of 
higher-order terms. Thus, the most important ingredient to the final stage of building the 
EFT is power counting. 

The controlled, natural expansion resulting from an EFT can be contrasted favorably 
with the characteristics of a theory constructed with traditional model-building strate- 
gies. A successful model takes a complex problem and reduces it reasonably to only a 
few degrees of freedom and types of interactions. This is a valid way to build a physical 
understanding of the problem, in that it allows the identification of the important sym- 
metries and degrees of freedom. It falls short in quantitative predictive power, however, 
because there is no way to place a bound on the error arising from effects not included in 
the model. EFT, on the other hand, provides a controlled expansion of the most general 
dynamics, so errors are well bounded and predictable. The value of such an advantage 
in comparing theoretical predictions to experiment cannot be overstated. 



A classical example 

At this point, some simple examples of EFT calculations will prove fruitful. We begin 
with a problem drawn from classical physics, a light object interacting via gravity with a 
much larger object, for example an apple near the surface of the earth. We can easily 
identify the relevant degree of freedom and symmetries of this simple problem. By 
experimenting with various objects close to the surface of the earth, we find that, to 
a good approximation, the degree of freedom is the mass m, while the symmetries 
involve translations parallel to the earth's surface and rotations about an axis normal 
to it. According to our recipe for an effective theory, we write down the most general 
potential reflecting these properties which is a power series in the height h: 

oo 

Veff{h)=mY,8ih'. (12) 

!=0 

Each of the g/s are parameters that can be fit to experimental data. (Here we are of course 
neglecting any quantum corrections, so there are no terms which are non-analytic in h 
and Vgyy is directly observable.) The first term is an irrelevant constant that depends on 
the arbitrary choice of the zero in energy. The second term is the familiar mgh, where 
g is the acceleration of a free body due to gravity as measured at h = 0. The linear 
form of the gravitational potential, however, is simply an approximation: higher-order 
terms are corrections that could in principle be extracted from careful measurements. 
This expansion will breakdown at some energy E^^^^^, which can be used to define a 
distance scale R through E , = mgR. This is a controlled expansion in h suitable for 



the low-energy regime E <^ E^^^^. The assumption of naturalness means 

where the cfs are coefficients of 0(1). This reduces to 

Thanks to Newton's consideration of large distances, in this case we know the un- 
derlying theory. A more accurate expression for the gravitational potential follows from 
Newton's law of universal gravitation, 

V(h) = -GMm^— , (15) 
^ ^ R + h 

where G is the gravitational constant, M is the mass of the earth, and R is its radius. For 
small h, we can expand Eq. (15) as 

V{h)^mgRYj,-iy-'l^-j , (16) 

where we have simplified the notation with the relation g = GM/R^. The underlying 
potential (16) indeed gives the effective gravitational potential (12) for h <^ R. Matching 
the two expressions, the coefficients of the effective potential are easily seen to be 

8i+i = (-l)'l- (IV) 

The naturalness assumption (14) is verified — all of the c/s are ±1 — with the added 
understanding of the scale R as the earth's radius. 

We have seen that the familiar linear form of gravitational potential energy can be 
treated as an effective theory of the more general theory given by Newton's law of 
gravitation. It is worth noting in passing, however, that Newton's law could also be 
viewed as an effective theory of a general relativistic formulation of gravitation whose 
effects are only visible at even larger energies. Thus, we can begin to develop a picture 
of nature as an "onion", with each successively smaller energy scale being described by 
an effective theory of the last. 



Non-relativistic QED 

While a simple example from classical mechanics is instructive, we are now ready to 
increase our insight into the nature of EFT by examining a quantum-mechanical problem 
that bring us closer to our nuclear-physics interest. Let us attempt to find the effective 
Lagrangian for QED at low energies. 



(further interactions) 




FIGURE 1. Schematic representation of a class of processes involving a low-momentum fermion. The 
fermion first interacts with a real or virtual photon and then propagates until it interacts further. 



Consider fermions described by a field \j/ of mass m and charge e, and assume the 
underlying theory to be given by the QED Lagrangian 

^QED = W[i{^-ie4) -m]Y+ ^F^yF^^ (18) 

where (^^v) is the photon field (field strength). 

What happens in processes where all particles have momenta of the same order 
Q <^ m? Consider the part of a diagram found in Fig. 1, where the fermion interacts 
with a low-energy photon, either real or virtual. In this non-relativistic regime we can 
write 

IpI ~ 1^1 ~ 2 < m, (19) 
/= y^2 + m2~m + o(^^j , (20) 

- 1^1 ~ Q. (21) 
The internal fermion line (with momentum p + q) contributes a factor to these diagrams 

4 + m) 



+....^l±f + (22) 



2pV + «e "' q^ + ie 2 

Here ". . ." stand for terms which are suppressed by powers of Q/m. 

The last equality in Eq. (22) contains the operator we will call P^, which can be 
defined together with the complementary P_, 

^± = -i^- (23) 



It is easy to verify that these operators possess the properties of idempotent projection 
operators: 

pI = P±, P±P^ = 0- (24) 

The presence of P_^_ in Eq. (22) indicates that the fermion propagation is given by 
the particle propagating forward in time. This suggests that we separate the nearly-inert 
mass m and the antiparticle component from the field, by rewriting the Xj/ ina "heavy- 
fermion formalism" [9] 

e"''"V± = P±V^. (25) 

The Lagrangian then becomes 



^QED = ¥+ido\j/+ - V-iy- + {if- ^V-)V+ + V-i-i^Q - 2m)v^_ + . . . (26) 

We calculate the low-energy Lagrangian from Eq. (10). Completing the square, doing 
the gaussian integral over X]/-, and renaming \]/_^_ to X]/ yields 

^eff (V^) = WK^o - ieAQ)xif+ ^ V7(V - /M) V+ . . . + ^F^yF^''' + ... (27) 

Of course, in the real world, there exist more than one particle coupling to photons, 
and other gauge bosons coupling to these particles. Because the photon is massless while 
weak bosons are not, the latter can be integrated out at low energies. Likewise, we can 
integrate out other, heavier fermions. This results in additional terms in the effective 
Lagrangian, such as a Pauli term that gives rise to an anomalous magnetic moment 
through a parameter /c: 

^f;t'^ (V') = i^^'iJkWO,WF,j + . . . (28) 

The new parameters can be calculated from the underlying theory. All interactions in 
£ (higher) suppressed by powers of a heavy mass. The Pauli term, for example, is 

oc 1 /m, and Jc is of 0(1) unless the particle represented by \j/ has only weak interactions. 

Now, according to Weinberg's "theorem", we do not need to go through this whole 
song and dance. If we directly construct the most general effective Lagrangian involving 
the fermion \j/ and photon A that is invariant under gauge transformations, parity, time 
reversal, and (non-relativistic) Lorentz boosts, we find the L^j^j above. 



Bound states in QED 

With the above effective Lagrangian, we can return to the problem of the hydrogen 

atom, considered in the Introduction. We have already seen that this problem possesses 
multiple energy scales, much like the nuclear-physics problem we wish to solve. Dis- 
cussing non-relativistic bound states in QED is a good warm-up for tackling nuclear 
bound states. 

For simplicity, let us consider the interaction of two non-relativistic fermions of 
the same type. We denote the initial (final) center-of-mass momentum p (p'), with 
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FIGURE 2. Diagrams representing all photon-exchange interactions between two heavy fermions in 
the center of momentum frame. The interaction T is shorthand for all interactions which can occur. Note 
that is a dummy 4-momentum, and as such should be integrated over in the two loop diagrams. 



\p\ ~ Ip'I ~ 2 and ~ p'^ ~ Q^/2m. Since heavy-fermion number is conserved, all 
diagrams that contribute to the scattering amplitude T have two fermion lines that go 
through. Let us consider first the diagrams made of photon exchange only, that is, of 
the type in Fig. 2, where the photon-fermion vertex comes from the first term in the 
Lagrangian (27). We want to estimate each to leading order in Q. 
The first-order diagram contributes a term 



{p - p'Y + ie (po - p'y -(p-p')2 + ie 
The second-order term in which the photon Unes cross is 
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The integral in can be evaluated as usual with a contour integral. It is clear that the first 

and second factors possess poles below the real axis, while the third and fourth terms 
each produce a pole above and below the real axis. We therefore minimize the algebra by 
closing in an infinite semicircle above the real axis: this way we avoid the two shallow 
poles from the nucleon propagators, and 
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(31) 



There is no surprise here: these contributions are, as ordinarily expected from loop 
diagrams, down by a factor of a compared to Eq. (29). 

Life is not so boring, fortunately. In the other second-order term, the photon lines do 
not cross. This contributes 
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The pole structure of Eq. (32) is the same as that of Eq. (30), except that the pole of the 
second factor lies in the upper half complex plane. In addition to residues from poles 
from photon propagators that contribute terms similar to the ones in Eq. (31), we cannot 
avoid a contribution from a shallow pole. 
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We note that 
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Thus, the term in r(^)(^^ from the shallow pole is enhanced over the other terms. Having 
done the time-component integrals in these diagrams, we can now think in terms of 
time-ordered perturbation theory. The reason time-ordered diagrams are useful here is 
that we are considering the interactions of non-relativistic particles. In time-ordered 
perturbation theory loops correspond to three-dimensional integrals. Each diagram in 
covariant perturbation theory unfolds into various time-ordered diagrams. The time- 
ordering of the four vertices in each second-order diagram is significant. Second-order 
terms, then, can be said to fall into two groups, those for which both photons are created 
and then both are destroyed, and those for which the first photon is destroyed before the 
second is created. These second diagrams are really iterations of first-order diagrams, 
and it is these that lead to the enhancement. It is clear that they can only occur as part of 
7^(2)(=)^ since r*^^'**^^) diagrams contain a crossing of photon lines by definition. 

Infrared enhancements will appear at all orders in perturbation theory, always from 
"two-fermion reducible" diagrams. The leading terms at each order form a series that 
goes roughly as 
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FIGURE 3. The Lippmann-Schwinger equation with potential V. 




FIGURE 4. Some time-ordered diagrams contributing to the two-particle potential V in non-relativistic 
QED. All orderings with at least one photon in intermediate states are included. A cross denotes energy 
corrections. 



It is clear from Eq. (35) that the series expansion breaks down at sufficiently small 
momenta such that m/Q compensates for a, that is, at 

Q ~ am, (36) 

indicating the existence of a bound state with binding energy 

q2. 

Br^^^ a^m. (37) 
m 

Hence we reproduce the estimates (5) and (6) arrived at earlier. 

This is no? just a more complicated way of generating known results: it is a more 
fundamental way. Of course, the actual resummation of this series can only be carried 
out by explicitly solving the Lippmann-Schwinger (or, equivalently, the Schrodinger) 
equation with a potential V , as represented in Fig. 3. The potential is given in leading 
order by the Coulomb potential given in Eq. (29). But from quantum mechanics alone 
we are at a loss about how to improve on the Coulomb potential. With non-relativistic 
QED, the prescription is clear, and it is illustrated in Fig. 4. The Coulomb potential of 
0{e^ /Q^) is the result from the lowest-order, static photon exchange. The corrections 
in the potential come from all other two-fermion irreducible time-ordered diagrams. 
As we have seem above, the irreducible second-order diagrams generate corrections of 
0{ae^ /Q^). The recoil part of the one-photon exchange (29) is an Q{Q/mx /Q^) 
correction, which for the bound state (where Q/m a) is comparable to two-photon 
exchange. And so on. 

A remarkable feature of the discussion above is that all the photon-exchange diagrams 
considered scaled with negative powers of Q. This is rooted on our temporary restriction 



to diagrams built from the simplest interaction in the Lagrangian (27). It means that no 
new ultraviolet cutoff is needed in loop integrals, once the bare fermion charge and mass 
are adjusted to remove cutoff dependence in diagrams involving a single fermion. At 
some point in the past great importance was attributed to this "renormalizability". Nowa- 
days, we realize that its simple form in QED is limited. Interactions with more deriva- 
tives, which are not forbidden by symmetry and thus exist, spoil this simplification. 
Piling up sufficiently many higher-derivative interactions from Eq. (28), a two-fermion 
diagram will scale as a positive power of Q, and diverge in the absence of a cutoff. How- 
ever, because four-fermion interactions can be constructed that are gauge invariant, the 
new cutoff dependence can still be eliminated by a shift in four-fermion parameters. The 
EFT is still renormalizable, in the sense that at any given power of Q a finite number 
of parameters will remove cutoff dependence from observables (renormalization-group 
invariance). QED can be deceivingly simple because the symmetries allow couplings of 
dimension four, which dominate. As we are going to see next, the same does not happen 
in the EFT of QCD. 



BUILDING A QCD-BASED EFT 
Standard model and QCD 

We now focus on the energy region of interest. The EFT corresponding to the standard 
model at an energy scale of a few GeV involves only the lightest leptons and quarks, 
gluons, and the photon; weak gauge bosons, and heavy leptons and quarks can be 
integrated out. For simplicity, here we focus on interactions involving the lightest u and 
d quarks. They can be arranged in a flavor doublet q = (^). The relevant pieces of the 
effective Lagrangian at this scale are 

3 1 

^ = -Y,ff{^-ie4 + mf)lf--F^yF^''' 
/=1 ^ 

-q{^ - igs^- ieQ4)q -:^{mu + m^)qq + :^{m^- mu)qT^q 
-^Tr[G^vG^'] + ^^e^.paTr[G^''GP-] + . . . (38) 

Here the /^^'s are the lepton fields with mass m^, (A^) is the gluon (photon) field of 
strength G^y i^iiv)^ is the usual Pauli matrix, Q= 1/6 + T3/2 is the quark charge 
matrix, and ". . ." denote higher-dimension terms. Higher-dimension interactions are 
suppressed by powers of the masses of the heavy particles that have been integrated 
out. We will neglect them in a (very good) first approximation. We will also neglect 
the theta term, since the strong CP parameter 9 is found to be unnaturally small (the 
so-called strong CP problem). The effect of these terms can be easily incorporated in a 
more comprehensive analysis. 

The leading Lagrangian (QCD + QED) becomes then invariant under parity and time- 
reversal transformations. The quark/gluon sector has 4 remaining parameters: the gauge 



couplings for strong (^.v) and electromagnetic (e) interactions, and the masses of the up 
(m„) and the down (m^) quarks. 

In the chiral limit, that is, when m„ = = and e = 0, the action of the QCD 

Lagrangian becomes invariant under scale transformations: A^/^^, and 

-^XGjx. Thus there is no dimensionful parameter in the Lagrangian; there is only one 
dimensionless parameter which is the strong-interaction coupling constant gv- However, 
the scale symmetry is explicitly broken due to quantum corrections that probe high 
energies. A regulator is necessary, thus introducing a dimensionful parameter in the 
problem. As a result, gs becomes a "running coupling" that is a function of a momentum 
scale. We observe in fact that gs grows as the scale is lowered, and reaches gs ~ 1 near 1 
GeV. Perturbation theory in gs, a useful tool at high energies, fails miserably at nuclear 
energies. 



Basic assumptions about a QCD-based EFT 

Thus in our pursuit to build an EFT, we have to take the properties of the nonper- 
turbative QCD into account, even though we cannot at present deduce them. The first 
assumption we make is that QCD "confines": only colorless states (hadrons) are asymp- 
totic states. In the low-energy EFT, only hadrons need to be incorporated as fields. Our 
second assumption is that QCD is natural. Since almost all hadron masses ^ 1 GeV, 
we conclude that there is a characteristic QCD mass scale Mg^^^ ~ 1 GeV. This sets the 
limit of validity of the EFT. 

Because we want to discuss nuclei, we certainly want to include the nucleon in the 
EFT. Because the delta-nucleon mass difference 5 = m^ — nij^ ~ 300 MeV ^ ^qcd^ 
if we want to explore the full region of energies accessible to EFT, we incorporate the 
delta isobar as well. Nucleon and delta fields form isospin multiplets as following: 



N={P,) and A= . (39) 




Below we denote these fermions by i//. Other baryon states such as A'^* are not included, 
as their mass differences to the nucleon approach M^^^ (m^* — ~ 500 MeV, etc.). 
Mesons such as the p and the (O are not included either because they have masses 
0{Mq^^). (Moreover, they interact with the pions and nucleons via dimension-four 
interactions that are not weak, which is at present an insurmountable obstacle for a 
systematic approach.) 

Does the pion mass, niji ~ 140 MeV <S ^qcd^ imply a breakdown of naturalness? 
No! We notice that in the chiral limit the leading-order Lagrangian (38) has an SUjil) x 
5C/^(2) ~ 50(4) chiral symmetry, since it is invariant under independent SU (2) rotations 
of left- and right-handed quarks, 

1r{l) = ^ ^Mi\^L) ■ ^) ^R{L) ■ (40) 



This symmetry is not manifest in the hadron spectrum. There is, for example, no scalar 
particle degenerate with the three pions. We thus make our third assumption, that chiral 
symmetry is spontaneously broken down to its diagonal subgroup, the SU^_^j^{2) ~ 
50(3) of isospin. Although we cannot at present calculate it, the effective potential of 
QCD, when plotted as function of quark bilinears {qxiy^q, qq), has to have a "Mexican- 
hat" shape. The symmetry of the potential is 50(4). The degenerate minima form thus 
a (four-dimensional) "chiral circle", the radius of which we call fj^. (Later, we find out 
that fji so defined coincides with the pion decay constant.) The symmetry is broken 
because the world actually sits on one particular point of this circle (which is selected 
by the quark-mass perturbation). Excitations orthogonal to the circle should have mass 
m<j ~ ^gcD' since this is the scale that sets the curvature of the effective potential. 
Excitations along the circle, on the other hand, are massless; the three Goldstone bosons 
can be associated with the light pions. 

Interactions of pions in the EFT have to reproduce these symmetry properties. There 
are several possible parametrizations of the chiral circle which correspond to different 
choices of pion fields K, such as the a-model-like, the Callan-Coleman-Wess-Zumino 
construction, and the stereographic projection. Since a chiral transformation rotates the 
chiral circle, the infinitesimal transformation of the pions is Tt ^ It + £. The EFT 
Lagrangian has to have the same symmetry, which implies that the pion fields can 
be chosen to couple derivatively. This is very important because it means that pion 
interactions involve the small momentum Q. Moreover, because the derivative is on the 
chiral circle, non-linear terms 

D-i = l-;rV4/^ + ... (41) 

appear. There is a well-defined procedure to construct the Lagrangian, the theory of non- 
linear realizations of a symmetry [10]. We first define co variant objects, such as covariant 
derivatives that involve the associated non-linear terms: the pion covariant derivative 
and the fermion covariant derivative D^. Then, the invariant Lagrangian can be 
constructed as the most general 50(3) -invariant Lagrangian made out of covariant 
objects. 

We can now relax the restriction to the chiral limit and consider the effect of non- 
vanishing quark masses. The important point is that the mass terms break 50(4) in a 
specific way. The common-mass term is the fourth component of an 50(4) vector, and 
the mass-difference term is the third component of another 50(4) vector. The common- 
mass term, for example, tilts the Mexican hat in the qq direction. The chiral circle is no 
longer degenerate, and the pions get a common mass m| ~ (m^ + m^^)Mg(-.^. In the EFT, 
we construct all operators that break chiral symmetry in the same way. Their coefficients 
will be proportional to powers of (m„-|-m^) and (m^ — m„) = e{mu + m^). 

Finally, electromagnetic interactions via "soft" photons can be added as well through 
all gauge-invariant terms. The one complication is that the integrating out of "hard" pho- 
tons generates operators that do not necessarily involve soft photons. Such "indirect" 
electromagnetic operators can be constructed by looking at the chiral transformation 
properties of the (non-local) four-quark operators they produce. One can show that these 
operators break chiral symmetry as 34 and 34-34 components of 50(4) antisymmet- 
ric tensors. The corresponding operators in the EFT will break isospin with strengths 



proportional to powers of e . 

To sum it up, the EFT Lagrangian has two classes of interactions. One class is chiral 
invariant; they involve powers of the momentum and no dimension-four interactions are 
permissible. Another class is chiral-symmetry breaking; dimension-four operators are 
among the ones allowed, but they are suppressed by powers of the small quark masses 
or charges. As a consequence, all pion interactions are weak at low energies. 



The chiral Lagrangian 

Based on the above, the most general Lagrangian takes the schematic form: 

f 



mnpqf 



xf^MhcD, (42) 



where the C^^p^j- are the unknown constants in the EFT. By naturalness these are of 

0(1) if isospin conserving and 0(e) and O(e^) if isospin breaking. 

It is convenient to introduce the "chiral index" A that counts the inverse powers of the 
large mass scale Mg^^^, 

f f 

A = n + m + q + 2p+^--2 = d+--2. (43) 

Because of chiral symmetry, there are no interactions with negative chiral index. We can 
then split the chiral Lagrangian in pieces L^^^: 

£(A)^ ^(0)^^(1) ^^(2) (44) 

A=0 

The form of the Lagrangian depends on the parametrization of the pion fields. For 
example, using stereographic coordinates, the lower-order Lagrangians are [11, 12, 13] 
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£(0) = -IflDl-'^D-^TC^+NUTiQN-lgjP-NhaN 



2 



N'^TSA + h.c. 



+ ..., (45) 
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iN'^TS-VA + h.c^+... 



(46) 



Here /i^ = 0(1) and B- = 0(l/A/g^^) are undetermined constants, to be obtained 
either by solving QCD or by fitting data; ". . . " stand for other terms involving the delta 
isobar. Higher-index interactions can be constructed similarly. 




Why is nuclear physics so interesting? 

Processes that involve at most one nucleon (A < 1) can be easily described in the EFT. 
If all external momenta are of the same order Q ~ M„„c, there are only two scales Q and 
^QCD- ^ generic contribution to an amplitude can be written as 

(47) 

where 5" is a dimensionless non-analytic function and X is a normalization factor. 
Counting powers of 2 in a particular diagram can be done as for the superficial degree 
of divergence. For a diagram with L loops and vertices of index A we find [7] 

V = 2-A + 2L + £V^A. (48) 

A 

This formula is important because chiral symmetry places a lower bound on the interac- 
tion index A > 0. Since L is bounded from below (L > 0), V > V^^-^ = 2 — A for strong 
interactions. An expansion in Q/Mq^^j^ results. It starts at v = V^-^ with tree (L = 0) 
diagrams built out of vertices of index (E^a^ ^ ^^^^ proceeds at v = V^^-^ -|- 1 
with further tree diagrams, now with one vertex of index 1, the remaining having index 
(L^A^ ~ These first two orders are equivalent to the current algebra of the 1960's, 
but now unitarity corrections can be accounted for systematically. At v = V^-^ + 2, for 
example, besides tree diagrams with one index-2 interaction or two index- 1 interac- 
tions (L V^-^ = 2), there are also one-loop (L — I) diagrams built out of index-0 vertices 
(L^A^ ~ This is generalized to higher orders in obvious fashion. An estimate of 
the expansion parameter is ~ nij^/nip ^ 0.2. In this context the EFT is called Chiral 
Perturbation Theory (ChPT). See Ref. [14] for a review. 



FIGURE 5. Simplest pion-exchange diagrams in the A^A'^ amplitude. 



Processes with two or more nucleons (A > 2) are much more interesting. One can see 
this by studying the simplest Feynman diagrams contributing to the A^A^ ampUtude T^^, 
shown in Fig. 5. The first diagram results in the term: 



(49) 



This is analogous to Coulomb-photon exchange in QED. By following the same steps 
as in that example, we see that the crossed-box diagram is 



fiUK)\0)\0)\0V\Q) AKflfV 



Jk 

which is small, given that 47cfj^ 
However, the box diagram is 



(50) 



QCD- 



The scale that sets the relative size of this diagram is not Mg^^ but 
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(51) 



M, 



NN 



(52) 



QCD- 



which numerically is ~ ~ Mnuc -C 

As in the QED case, the infrared enhancement of m^y / Q over irreducible states appears 
in all reducible intermediate states. Contrary to ordinary ChPT, we need to resum the 
leading terms; in a schematic way. 
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(53) 



Nuclear bound states thus appear at 

e-M^^, (54) 

with binding energy 



2 Ml,, 

B^^ ^. (55) 



NN 



Numerically this is ~ Mq^j^/ [AnY a; 10 MeV, which explains why nuclei are so shallow 
when compared to the characteristic QCD scale. 



Deriving a potential for nuclear physics 

We can now carry out systematic calculations by finding a power counting for the 
potential, defined as the sum of all irreducible diagrams, and then solving the Lippman- 
Schwinger equation (Fig. 3) order by order. 

Weinberg [15] was the first to suggest a simple power counting. He reasoned that, 
because the enhanced states are by construction removed from the potential, powers of 
Q could be counted in the potential as in ChPT. This results in an expansion of the form: 

oo 

V = £ c^Q\ (56) 
v=0 

where 

v = 4-A + 2(L-C)+£V^A. (57) 

A 

The only difference compared to Eq. (48) is the number of separately connected pieces 
C > 1. They arise because the potential is just part of the full amplitude and does not 
need to be fully connected. 

According to this power counting, in leading order we have, besides one-pion ex- 
change (OPE) between two nucleons, also non-derivative two-nucleon contact interac- 
tions (the C^^^ and C^^^ terms in Eq. (45).) A calculation of all (isospin-conserving) 
contributions to the two-nucleon potential up to V = v^^^ -\- 3 was carried out in Refs. 
[11, 12] using time-ordered perturbation theory. Some diagrams are shown in Fig. 6. 
To this order, the potential has all the spin-isospin structure of the phenomenological 
models, but its profile is determined by explicit degrees of freedom, symmetries, and 
power counting. The power counting suggests a hierarchy of short-range effects: S waves 
should depend strongly on the short-range parameters C^'^^ ; contact interactions affect 
P-wave phase shifts only in subleading order, so their effect should be smaller and ap- 
proximately linear; D waves are affected by contact interactions only via mixing, while 
higher waves should be essentially determined by pion exchange. Chiral symmetry is 
particularly influential in the two-pion exchange (TPE) piece. The latter includes a par- 
ticular form of terms previously considered [16], plus a few new terms. Those terms 
involving the 5-'s and the deltas provide the only form of correlated TPE to this order, 
as graphs where pions interact in flight appear only in next order and should thus be 



t: 



FIGURE 6. Some time-ordered diagrams contributing to the two-nucleon potential V^^ in the EFT. 
(Double) solid lines represent nucleons (and/or deltas), dashed lines pions, a heavy dot an interaction in 
a dot within a circle an interaction in iL''', and a dot within two circles an interaction in . First 
line corresponds to V = V^;„, second and third lines to V = V^;„ + 2, fourth line to v = V^;„ + 3, and ". . ." 
denote v > v^,,,,, + 4. AU orderings with at least one pion or delta in intermediate states are included. Not 
shown are diagrams contributing only to renormalization of parameters. 



relatively small. The sum of the term and the corresponding delta term (related to the 
nucleon axial polarizability), is particularly important in providing an isoscalar central 
force. Not surprisingly, in the chiral limit these potentials behave at large separations as 
van der Waals forces. Isospin- violating pieces of the potential can be calculated as well 
[19]. 

Using this potential, one can then solve the Schrodinger equation. For each cutoff 
the parameters are fitted to data at selected energies. A sample of the results [12] for the 
lowest, most important partial waves is presented in Fig. 7, together with phases from the 
Nijmegen phase-shift analysis [17]. Quality of the fits is typical of other waves. Waves 
higher than F-waves were found to be mostly described well by pion exchange alone, as 
expected. Deuteron quantities are shown in Table 1. Electromagnetic quantities refer to 
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FIGURE 7. ^Sq, and ^D^ NN phase shifts and mixing angle in degrees as functions of the 
laboratory energy in MeV: EFT up to v = 3 for cutoffs of 500 (dotted), 780 (dashed), and 1000 MeV 
(solid line); and Nijmegen PSA (squares). 



the contributions from lowest-order yA^A^ couplings only, not to a consistent calculation 
which would include sub-leading one- and two-nucleon effects. The predicted 5-wave 

scattering lengths (not used to constrain the fit) were found to be a^^^^ — —15.0 fm 

and ~ 5.46 fm. Important central attraction comes from the B-'s and deltas, and 
indeed the central potential does resemble those from models that include o and (O 
meson exchange explicitly [18] (see Fig. 8). Thus the properties of the mesons extracted 
from phenomenological models have limited meaning. Values for the parameters are 
listed in Ref. [12]. Reasonable values were found for quantities known at the time, for 
example = 1.33 (in agreement with the Goldberger-Treiman relation) and /j^ = 2.03 
(smaller but not too far off the large-X^ value). However, the values for the 5-'s came 
out somewhat different than those found later in nN scattering. For a cutoff of A = 780 
MeV, the coefficients of the contact interactions were found to scale approximately 
with M^^n~ 500 MeV. 
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FIGURE 8. The isoscalar central potential generated by two-pion exchange compared with phenomeno- 
logical (7 + 0) contributions [18]. 



The EFT also offers some insight into few-nucleon forces. Weinberg's power counting 
embodied in Eq. (57) suggests a hierarchy of few-nucleon forces. In leading order 
(v = V^^^ = 6 — 3A), C is maximum, so we have only pairwise interactions via the 
leading two-nucleon potential. We can easily verify that, if the delta is kept as an 
explicit degree of freedom, a 3N potential will arise at V = V^-,^ + 2, a 4N potential 
at V = V^-^ + 4, and so on. It is (approximate) chiral symmetry therefore that implies 
that n-nucleon forces V^^y obey a hierarchy of the type 



with (y„^) denoting the contribution per n-plet. If 1(^2^)1 ~ MQfjjy/{4n)^ a; 10 MeV, 



TABLE 1. Results from EFT fits at v = 3 for various 
cut-offs A and experimental values for the deuteron bind- 
ing energy (B), magnetic dipole moment (;U^), electric 
quadrupole moment (Q^), asymptotic D/S ratio (tj), and 
D-state probability (f^). 

Deuteron A (MeV) 

quantities 500 780 1000 Experiment 

B(MeV) 2.15 2.24 2.18 2.224579(9) 

Aij(jU^) 0.863 0.863 0.866 0.857406(1) 

e£(fm2) 0.246 0.249 0.237 0.2859(3) 

ri 0.0229 0.0244 0.0230 0.0271(4) 
Pr,(%) 2.98 2.86 2.40 
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FIGURE 9. Some time-ordered diagrams contributing to the 3A^ potential in the pionful EFT. (Double) 
solid lines represent nucleons (and/or deltas), dashed lines pions, a heavy dot an interaction in iL'"', and a 
dot within a circle an interaction in L^^' . First line corresponds to V = V^,„ + 2, second line to V = V^,„ + 3, 
and ". . ." denote v > v^;„ + 4. All nucleon permutations and orderings with at least one pion or delta in 
intermediate states are included. 



we can expect | (^3^) | ~ 0.5 MeV, | {V^^^) \ ~ 0.02 MeV, and so on. This is in accord with 
detailed few-nucleon phenomenology based on potentials that include small 3^" and no 
AN forces. This is shown in Table 2 in the case of the AV18/IL2 potential [20]. 

The new forces that appear in systems with more than two nucleons have been derived 
in Refs. [11, 13]. The relevant terms up to + 3 are shown in Fig. 9. The leading 
3A^ potential has components with three different ranges: TPE; OPE/short; and purely 
short. The TPE part of the potential is determined in terms of kN scattering observables, 
and is similar to existing two-pion-exchange 3N potentials [21]. The novel OPE/short- 
range components involve 7r{N^N)^ interactions of strengths that are not fixed by chiral 
symmetry alone but can be determined from reactions involving only two nucleons, 
such as A^A^ — > A^A^;r [22]. The purely short-range components of the potential can only 
be determined from few-nucleon systems. 

The nuclear potential from EFT has been further elaborated in the last couple of years. 
Calculations are being pushed to next order [23], better fits (with more extensive input 
from ^rA'^ scattering) to A^A'^ data have been achieved [24], and the first results for 3A'^ and 



TABLE 2. Contributions of the two-, three- and four-nucleon potentials (per doublet, per triplet, 
and per quadruplet, respectively): Weinberg power counting (W pc) and calculations with the 
AV18/IL2 potential for the ground states of various light nuclei (^H, ^H, etc.). 
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FIGURE 10. The short-ranged (p, etc.) cloud around the nucleon and the long-ranged {%) cloud 
surrounding it. 



4A'^ systems have appeared [25]. These developments are all reviewed in Ref. [1]. 

Weinberg's power counting thus provides a rationale to understand much of the 
phenomenology of light nuclei. At low energies, the nucleon can be visualized as in 
Fig. 10: the light pion states form a sparse outer cloud, leading to a loop expansion in 
Q/AkFjj:, while the high-energy states form a dense inner cloud, amenable to a multipole 
expansion in terms of Q/nip. Since most of the time there is only a single pion "in the 
air", in the nucleus the interaction among nucleons is mostly pairwise, resulting also in 
a cluster expansion for the potential. 



Renormalization 

Despite all similarities, there is an important difference between the nuclear and 
atomic potentials. Because of chiral symmetry, the potential from one-pion exchange 
goes as Q^, from two-pion exchange as Q^, and so on. The loops generated by the 
iteration of these potentials are then sensitive to the cutoff. In coordinate- space parlance, 
these potentials are singular, that is, they behave for small radial distances r as 1/r^, 
and so on. The Schrodinger equation for such potentials is not well defined. 
Clearly, one cannot consider pion exchange without including enough short-distance 
operators to erase the cutoff dependence from observables. 

Now, Weinberg's power counting automatically predicts which contact operators have 
to be taken at each order. Following this prescription, it does seem that the fits are of the 
same quality for different cutoffs. However, the numerical character of the calculations 
makes it hard to offer a proof of consistency. It is not apparent, for example, how a 
momentum-independent short-range interaction can renormalize the OPE potential (49). 



Much effort has been devoted to this issue in the last five years [1]. A search for more 
analytical approaches led to extensive studies of the EFT at momenta Q <^ ni-jz- In this 
regime, pions can be treated as heavy particles. All interactions are of contact type, and 
form an expansion in Q/Mnuc- This "pionless" EFT can be solved by hand in the A^A^ 
sector, and by very simple numerical means in the 3N and 4A'^ systems. Power counting 
and renormalization are not without surprises, but they can be done consistently. Of 
course, the pionless EFT does not address the problem of pion exchange directly, but it 
goes a long way in elucidating how EFTs work in a non-perturbative context. Moreover, 
it is useful for many low-energy reactions, to which it has by now been applied. 

As Q is increased, the approximation of zero-range propagation for pions becomes 
less and less reliable. At momenta Q ~ m^^, pions have to be included explicitly in the 
theory. Once pions have been reinstated, we can still consider the low-Q region. Because 
of chiral symmetry, it is reasonable to suppose that pion interactions are perturbative at 
sufficiently low Q. Based on the power counting of the pionless EFT, a new power 
counting was formulated that led to a manifestly consistent EFT where pion exchange 
was treated in perturbation theory [26]. Unfortunately, it has been shown that the range 
of validity of this power counting is Q ^ 100 MeV [27]. In this region the simplest 
pionless EFT is more useful. 

Recently, the non-perturbative renormalization of pion exchange was analyzed [28]. 
It was found that OPE can be renormalized a la Weinberg, provided one expands in the 
pion mass even the pion propagators. If one does not do expand, one makes numerically 
small but conceptually large errors. The full implication of this conclusion to previous 
results remains to be explored. 



Nuclear matter on a lattice 

The nuclear EFT that has been sketched above can be used to study nuclear matter 
at finite temperature. This provides a link with RHIC physics that is the subject of this 

meeting. 

The general concept of a nuclear matter calculation consists of nucleons interacting 
via a variety of components of the nuclear potential. While the ultimate goal is to use 
EFT interactions, let us concentrate for simplicity on few parts of the A^A^ potential, 
namely central, spin- and isospin-exchange. The Hamiltonian, 

f( = X + V, (59) 

can be expressed in second quantization and contains kinetic and potential operators. 
The kinetic term is written as 

±=-Va1--^V', (60) 
2m^ 

while the potential is taken as 

V = Co(v.V)' + C2(v^Vv^)'. (61) 



The fermion operator i//^ creates a nucleon of spin and isospin (c7, t) at location x, while 
its adjoint y/ destroys it. This Hamiltonian describes merely a toy model for developing 
the formalism. More complicated Hamiltonians arising from full-fledged nuclear EFT 
can be later used for more realistic calculations. 

How can we deal with the complication of many nucleons? A natural approach, 
patterned after similar attempts in QCD, is to investigate nuclear matter on a three- 
dimensional cubic lattice of spacing a and periodic boundary conditions. We describe 
the nuclear-matter Monte Carlo method [29], which consists of the thermal formalism to 
express the grand-canonical partition function as an integral over single-body evolution 
operators. At its center stands the Hubbard-Stratonovitch transformation, which is used 
to reduce the many-body problem to an effective one-body problem. 

In order to study thermal properties of nuclear matter, the grand-canonical partition 
function at a given temperature T = j8~^ needs to be determined: 



Z = Tr 



exp (-j8 vaVv^))] = Tr [U] , (62) 



with jU as the chemical potential. U is called the imaginary-time evolution operator of the 
system and is a many-body operator; the trace is taken over all many-body states. The 
partition function Z is an exponential over all one- and two-body operators (and therefore 
interactions) present in the system. It is impossible to deal with Z in this form, because 
the number of many-body correlations that have to be kept track of grows rapidly with 
system size. We therefore find an expression for Z that is based on a single-particle 
representation, replacing the many-body problem with that of non-interacting nucleons 
that are coupled to a heat bath of auxiliary fields. A thermal observable (O) is then 
expressed as [30] 



where the Gaussian factor G is given by 

^(^) = ft nn^'^p {-\^i\xhn,) ■ (64) 

m=l Xn i 

The integrals in this equation can then be evaluated using the Metropolis algorithm [31] 
(Monte Carlo simulations). 

Fig. 11 [29] shows that such calculations are feasible. It displays the energy per 
particle E/Afor symmetric nuclear matter and for pure neutron matter as a function of 
density p and for different temperatures. The two interaction parameters were adjusted 
so that the known qualitative features of the symmetric nuclear matter are reproduced. 
With decreasing temperature, symmetric nuclear matter develops a minimum at p = 
0.32 fm~^ first, which is most pronounced between 10 — 14 MeV, before it shifts to 
lower densities. At T = 3.3 MeV and T = 5.9 MeV the minimum is very broad, making 
matter softer. For high temperatures and/or high density, the simulation suffers from 
the fact that it runs out of model space. At T = 50 MeV the system behaves almost 
like a Fermi gas and the energy per particle should behave like ~ p^/^. Yet, the curve 




bends down. Also, for all other temperatures, the curves converge to the energy of the 
full lattice state, E /A — 5.96 MeV, as density increases. For sub-saturation densities 
the model gives more binding if compared to other calculations (see, for example, Ref. 
[32]), and the energy is not as high for densities beyond saturation. At p = 0.32 fm^^, 
E /A as a function of temperature has a minimum at T ?s 10 MeV which means that 
at even lower temperatures E/A increases again. We see some evidence for a liquid- 
gas phase transition. The uncertainties for pure neutron matter are much larger than for 
symmetrical nuclear matter. As a potential, we used the parameters obtained from the fit 
to symmetric nuclear matter, even though we could have fitted the potential parameters 
for this case anew. Therefore we view the results for pure neutron matter more as a 
test to see how well the given potential already reproduces the energy. Note that the 
slopes of the curves at high temperatures are not negative as they are for symmetric 
nuclear matter. But clearly, we cannot conclude that the energies at T = 3.3 MeV have 
converged to that of the ground state because the curve differs quite a bit from that of 
r = 5.9 MeV. At the lowest temperature they are 4 — 5 MeV higher than those of the 
ground state as calculated in Ref. [32], but the general shape of the curve is very similar. 
This is no surprise, since pure neutron matter is like a Fermi gas, with attractive forces 
between neutrons lowering the energies with respect to the non-interacting system. The 
search for any kind of phase transition in the range of 5 — 50 MeV was to no avail. 



CONCLUSION 

We have presented an introduction to EFT and some of the mileposts along the sinuous 
road from QCD to nuclear physics. We discussed some of the essential features of the 
interactions among nucleons and some of their consequences to few-nucleon systems. 
We can already see how to approach the solution of the important problem of finite- 
temperature nuclear matter. Yet, there is no denying that much is still to be done in 
straightening the path already traveled, and in breaking new ground. It is a good time to 
be on the road! 



ACKNOWLEDGMENTS 



The authors acknowledge the support and dedication of all PASI sponsors and organiz- 
ers, in particular Jan Rafelski, Bob Thews, Thomas Elze, Erasmo Ferreira, and Takeshi 
Kodama. UvK thanks the Nuclear Theory Group at the University of Washington for 
hospitality while part of this work was carried out, and RIKEN, Brookhaven National 
Laboratory and to the U.S. Department of Energy [DE-AC02-98CH 10886] for provid- 
ing the facilities essential for the completion of this work. This research was supported 
in part by a DOE Outstanding Junior Investigator Award (UvK). 



REFERENCES 

1. Bedaque, P.F., and van Kolck, U., nucl-th/0203055, to appear in Ann. Rev. Nucl. Part. Sci.; van 
Kolck, U., Nucl. Phys. A699, 33 (2002). 

2. Phillips, D., nucl-th/0203040, to appear in Czech. J. Phys. 

3. Gerber, P., and Leutwyler, H., Nucl. Phys. B321, 387 (1989). 

4. Hagedorn, R., and Rafelski, J., Phys. Lett. B97, 136 (1980). 

5. Karsch, E, AIP Conf. Proc. 602, 323 (2001). 

6. Braaten, E., and Nieto, A., Eur. Phys. J. Bll, 143 (1999); Hammer, H.-W., and Eumstahl, R.J., Nucl. 
Phys. A678, 277 (2000). 

7. Weinberg, S., Physica A96, 327 (1979). 

8. Ball, R.D., and Thorne, R.S., Ann. Phys. 236, 117 (1994). 

9. Georgi, H., Phys. Lett. B240, 447 (1990). 

10. Coleman, S., Wess, J., and Zumino, B., Phys. Rev. 177, 2239 (1969); Callan, C.G., Coleman, S., 

Wess, J., and Zumino, B., Phys. Rev. Ill, 2247 (1969). 

11. Ordonez, C, and van Kolck, U., Phys. Lett. B291, 459 (1992). 

12. Ordonez, C, Ray, L., and van Kolck, U., Phys. Rev. C53, 2086 (1996). 

13. U. van Kolck, Phys. Rev C49, 2932 (1994). 

14. Bernard, V, Kaiser, N, and MeiBner, U-G., Int. J. Mod. Phys. E4, 193 (1995). 

15. Weinberg, S., Phys. Lett. B251, 288 (1990); Nucl Phys. B363, 3 (1991). 

16. Brueckner, K., and Watson, K., Phys. Rev 92, 1023 (1953); Sugawara, H., and Okubo, S., Phys. Rev 
111, 605 (1960); Sugawara, H., and von Hippel, R, Phys. Rev 112, 1764 (1968). 

17. Stoks, V.G.J. , Klomp, R.A.M., Rentmeester, M.C.M., and de Swart, J.J., Phys. Rev C48, 792 (1993). 

18. Kaiser, N., Brockmann, R., and Weise, W., Nucl. Phys. A625, 758 (1997). 

19. van Kolck, U., Few-Body Syst. Suppl. 9 444 (1995); van Kolck, U., Friar, J.L., and Goldman, T., 
Phys. Lett. B371, 169 (1996); van Kolck, U., Rentmeester, M.C.M., Eriar, J.L., Goldman, T., and de 
Swart, J.J., Phys. Rev Lett. 80, 4386 (1998); Eriar, J.L., and van Kolck, U., Phys. Rev C60, 034006 
(1999); Niskanen, J.A., Phys. Rev C65, 037001 (2002). 

20. Pieper, S.C., and Wiringa R.B., Ann. Rev Nucl. Part. Sci. 51, 53 (2001). 

21. Friar, J.L., Huber, D., and van Kolck, U., Phys. Rev C59, 53 (1999). 

22. Huber, D., Friar, J.L., Nogga, A., Witala, H., and van Kolck, U., Few-Body Syst. 30, 95 (2001); 
Hanhart, C, van Kolck, U., and Miller, G.A., Phys. Rev. Lett. 85, 2905 (2000). 

23. Kaiser, N., Phys. Rev C61, 014003 (2000); C62, 024001 (2000); C63, 044010 (2001); C64, 057001 
(2001); C65, 017001 (2002). 

24. Epelbaum, E., Glockle, W., and MeiBner, U.-G., Nucl. Phys. A671, 295 (2000); Entem, D.R., and 
Machleidt, R., Phys. Lett. B524, 93 (2002); Walzl, M., MeiBner, U.-G., and Epelbaum, E., Nucl. 
Phys. A693, 663 (2001). 

25. Epelbaum, E., Kamada, H., Nogga, A., Witala, H., Glockle, W., and MeiBner, U.-G., Phys. Rev. Lett. 
86, 4787 (2001). 

26. Kaplan, D.B., Savage, M.J., and Wise, M.B., Nucl. Phys. B534, 329 (1998). 

27. Fleming, S., Mehen, T., and Stewart, I.W., Nucl. Phys. Mil, 313 (2000). 

28. Beane, S.R., Bedaque, RE, Savage, M.J., and van Kolck, U., Nucl. Phys. A700, 377 (2002). 



29. Miiller, H.-M., Koonin, S.E., Seki, R., and van Kolck, U., Phys. Rev. C61, 044320 (2000). 

30. Loh Jr., E.Y., and Gubernatis, J., in Electronic Phase Transitions, ed. by W. Hanke and Yu.V. Kopaev 
Elsevier, New York, 1992; Lang, G., Johnson, C, Koonin, S., and Ormand, W., Phys. Rev. C48, 1518 
(1993). 

31. Metropolis, N., Rosenbluth, A., Rosenbluth, M., Teller, A., and Teller, E., /. Chem. Phys. 21, 1087 
(1953). 

32. Wiringa, R.B., Fiks, V., and Fabrocini, A., Phys. Rev. C38, 1010 (1988); Akmal, A., and Pandhari- 
pande, V.R., Phys. Rev. C56, 2261 (1997). 



